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Abstract. We show that in the equatorial plane of marginally stable thick discs (with uniformly distributed specific angular 
momentum £(r, 6) = const) the orbital velocity relative to the LNRF has a positive radial gradient in the vicinity of black holes 
with a > 0.99979. The change of sign of the velocity gradient occurs just above the center of the thick toroidal discs, in the 
region where stable circular geodesies of the Kerr spacetime are allowed. The global character of the phenomenon is given 
in terms of topology changes of the von Zeipel surfaces (equivalent to equivelocity surfaces in the tori with £(r, 0) = const). 
Toroidal von Zeipel surfaces exist around the circle corresponding to the minimum of the equatorial LNRF velocity profile, 
indicating a possibility of development of some vertical instabilities in those parts of marginally stable tori with positive 
gradient of the LNRF velocity. Eventual oscillatory frequencies connected with the phenomenon are given in a coordinate- 
independent form. 
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1. INTRODUCTION 

High frequency (kHz) quasi-periodic oscillations (QPOs) with frequency ratios 3:2 (and sometimes 3:1) are observed 
in microquasars (see, e.g., [1, 2]). The same frequency ratios of QPOs in mHz are observed in the Galactic Center 
black hole Sgr A* [3, 4]. It is assumed now that the QPOs are related to the parametric or forced resonance [5] of the 
radial and vertical epicyclic oscillations or one of the epicyclic and orbital oscillations in accretion discs [6, 7, 8]. 

The oscillations could be related to both the thin Keplerian discs [9, 10] or the thick, toroidal accretion discs 
[11, 12]. The parametric resonance of the radial and vertical oscillations in the thin discs can explain the QPOs with the 
cog /co, - = 3 : 2 frequency ratio observed in all the microquasars and can put strong limits on the rotational parameter 
of their central black holes related to the limits on their mass [13]. 

Aschenbach [4] conjectured relation between the 3:2 and 3:1 resonance orbits by relating their Keplerian orbital 
velocities at an d r^ : i to be Q.^(r^ \;a) = 3Q.^(ry2',a), fixing thus the rotational parameter of black holes at the 
value of a = af = 0.99616. Further, he proposed that excitation of the oscillations at r = r^i can be related to two 
changes of sign of the radial gradient of the Keplerian orbital velocity as measured in the locally non-rotating frame 
(LNRF) [18] that occurs in vicinity of r = ry i for black holes with a > 0.9953. While the assumption of frequency 
commensurability of Keplerian orbits at rj, i and seems to be rather artificial because distant parts of the Keplerian 
disc have to be related, we consider the positive radial gradient of orbital velocity in LNRF nearby the ry \ orbit around 
black holes with a > 0.9953 to be a physically interesting phenomenon, even if a direct mechanism relating this to 
triggering of the excitation of radial (and vertical) epicyclic oscillations is unknown. We also show (Table 1) that 
physically relevant, i.e. coordinate-independent, frequencies, given by maximal positive radial gradient of the LNRF- 
orbital velocity, take locally values from tens to hundreds of Hz and for stationary observers at infinity tens of Hz (for 
typical stellar-mass black holes M ~ 10M Q ), which likely disables possibility to identify these frequencies directly 
with those of high frequency QPOs [13]. 

Because the accretion-disc regime will vary from thin Keplerian disc to thick toroidal disc with variations of 
accretion flow, we shall study here, without addressing details of the mechanism, whether the orbital velocity in 
LNRF can have positive gradient also for matter orbiting black holes in marginally stable thick discs with uniform 
distribution of the specific angular momentum (£(r, 9) = const), leading to a possibility to excite oscillations in the 
thick-disc accretion regime. Note that the assumption of uniform distribution of the specific angular momentum can 
be relevant at least at the inner parts of the thick disc and that matter in the disc follows nearly geodesic circular orbits 
nearby the centre of the disc and in the vicinity of its inner edge determined by the cusp of its critical equipotential 
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surface, see [14, 15]. 

In thick tori, it is necessary to have information about the character of the Aschenbach's phenomenon also outside the 
equatorial plane. We shall obtain such information by introducing the notion of von Zeipel radius analogical to the 
radius of gyration p introduced for the case of Kerr spacetimes in the framework of optical geometry by Abramowicz 
et al. [16], generalizing in one special way the definition used for static spacetimes [17]. The von Zeipel radius is 
defined in such a way that for the marginally stable tori the von Zeipel surfaces, i.e., the surfaces of constant values of 
3%, coincide with surfaces of constant orbital velocity relative to the LNRF. 



2. TOROIDAL MARGINALLY STABLE ACCRETION DISCS 

The perfect fluid stationary and axisymmetric toroidal discs are characterized by 4-velocity field U^ 1 — (£/', 0, 0, U' p ) 
with U' = U*{r, 9), = U^{r, 9), and by the distribution of specific angular momentum I = —Uy/Ut. The angular 
velocity of orbiting matter, Q. = U^/U', is then related to I by the formula 

£l = - l S»+8t 9 (1) 
£gt<p+g<pcp 

The marginally stable toroids are characterized by the uniform distribution of specific angular momentum £ = 
£(r, 9) = const and they are fully determined by the spacetime structure through equipotential surfaces of the potential 
W = W(r, 9) defined by the relation [14] 

w-w ra =in-p-, ( Ut f= ; (2) 

{Ut )in gtt£ Z + 2gt<pe + g<p<p 

the subscript "in" refers to the inner edge of the disc. 

In the Kerr spacetimes with the rotational parameter assumed to be a > 0, the relevant metric coefficients in the 
standard Boyer-Lindquist coordinates read: 

A -a 2 sin 2 9 2arsin 2 9 Asin 2 

gtt = ^ , gt<p = ^ , g<p<p = ^ , (3) 



where 

A = r 2 -2r + fl 2 , £ = r 2 + a 2 cos 2 9, A = (r 2 +a 2 ) 2 - Aa 2 sin 2 9. (4) 

The geometrical units, c = G = 1, together with putting the mass of the black hole equal to one, M = 1, are used to 
obtain completely dimensionless formulae hereafter. The relation (1) for the angular velocity of matter orbiting the 
black hole acquires the form 

-2 ■„ 2 a\e i o„„„;„2 , 



£l = £l(r,9;a,e) = ± „/ N _.„, 2 _ (5) 



(A-a 2 sm z 9)e + 2arsm z 9 
(A - liar) sin 2 9 

and the potential W, defined in Eq. (2), has the explicit form 



W = W(r, 9;a,£) = -In— , EA f" ® = — -. (6) 

V y 2 (r 2 + fl 2 -fl^) 2 sin 2 0-A(£-flsin 2 0) 2 

3. THE ORBITAL VELOCITY IN LNRF 

The locally non-rotating frames are given by the tetrad of 1 -forms [18] 

e W=^y /2 d* ) eW=(|y 2 dr, e( e )=£ 1 / 2 d0 ! e« = (^j^ S m9(d<p - codt), (7) 

where 

& = _gu L = 2ar 
A 
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is the angular velocity of LNRF. The azimuthal component of 3 -velocity in LNRF reads 



LNRF uV( y <»)■ W 

Substituting for the angular velocities i2 and © from the relations (5) and (8), respectively, we arrive at the formula 



_ A(A-a 2 sin 2 0)+4fl 2 r 2 sin 2 
!Ly/A(A-2a£r) sin0 



We focus our investigation to the motion in the equatorial plane, = n/2. Formally, this velocity vanishes for r — > °° 
and r — > r + = 1 + vl — a 2 , where the event horizon is located, i.e., there must be a change of its radial gradient for 
any case of values of the parameters a and I, contrary to the case of Keplerian orbits characterized by the Keplerian 
distributions of the angular velocity and the specific angular momentum 

1 r 2 ~2ar x l 2 + a 2 

fl = fl K (r; a )= e = tK(r,a)= f , (11) 

where the azimuthal component of the 3-velocity in LNRF in the equatorial plane, 9 = n/2, reads 

r (<P) (rj) (r 2 +a 2 ) 2 -a 2 A-2ar(ry 2 + a) 

and formally diverges at r — r+. 

Of course, for both thick tori and Keplerian discs we must consider the limit on the disc extension given by the 
innermost stable orbit. For Keplerian discs this is the marginally stable geodetical orbit, while for the thick tori 
this is an unstable circular geodesic keeped stable by pressure gradients and located between the marginally bound 
and the marginally stable geodetical orbits, with the radius being determined by the specific angular momentum 
£ = const € (/ ms ,/ m b) through the equation I = £jt(r; a); £ ms (4nb) denotes specific angular momentum of the circular 
marginally stable (marginally bound) geodesic. 

The radial gradient of the equatorial orbital velocity of thick discs reads 

= [A+(r-l)r}[r(r 2 + a 2 )-2a(e-a)}-r(3r 2 + a 2 )A £ 
dr ~ [r{r 2 + a 2 )-2a(£-a)] 2 VA 

so that it changes its orientation at radii determined for a given I <G (/msj^mb) by the condition 

9 ( \- i r 2 \{r 2 +a 2 ){ r -\)-2rA] 

£ = £ eK (r;a) = a + 2a[A + r(r _ 1}] • (14) 

Detailed discussion [19] shows that two changes of sign of d'f^'/dr can occur for Kerr black holes with the 
rotational parameter a > a c ( t hi c k) = 0.99979. The interval of relevant values of the specific angular momentum 
l € (^ms(fl))^ex(max)( a )) grows with a growing up to the critical value of a c ( m b) = 0.99998. For a > a c t m ^\, the interval 
of relevant values of £ € (£ m s(a) , £mb(a)) is narrowing with growing of the rotational parameter up to a = 1, which 
corresponds to a singular case where £ ms (a = 1) = £ m b{ a = 1) = 2. Notice that the situation becomes to be singular 
only in terms of the specific angular momentum; it is shown (see [18]) that for a — 1 both the total energy E and 
the axial angular momentum L differ at r ms and r m b, respectively, but their combination, £ = L/E, giving the specific 
angular momentum, coincides at these radii. 



4. VON ZEIPEL SURFACES 

It is useful to find global characteristics of the phenomenon that is manifested in the equatorial plane as the existence of 
a small region with positive gradient of the LNRF velocity. A physically reasonable way of defining a global quantity 
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FIGURE 1. Von Zeipel surfaces (meridional sections). For a < a c (bh) and any I, only one surface with a cusp in the equatorial 
plane and no closed (toroidal) surfaces exist. The cusp is, however, located outside the toroidal equilibrium configurations of perfect 
fluid. For a > a c (t,h) and ^ appropriately chosen, two surfaces with a cusp, or one surface with both the cusps, together with closed 
(toroidal) surfaces, exist located always inside the ergosphere (dashed surface) of a given spacetime. Moreover, if a > a c ( t hi c k)' 
both the outer cusp and the central ring of closed surfaces are located inside the toroidal equilibrium configurations corresponding 
to marginally stable thick discs (light-gray region; its shape is determined by the critical self-crossing equipotential surface of 
the potential W given by (6)). The cross (+) denotes the center of the torus. Dark region corresponds to the black hole. Figures 
illustrating all possible configurations of the von Zeipel surfaces are presented in [19]. Here we present the figure plotted for 
the parameters a = 0.99998, I = 2.0065. Critical value of the von Zeipel radius corresponding to the inner and the outer self- 
crossing surface is ^ c (; n ) = 3.429 and ,^ c ( out ) = 3.804, respectively, the central ring of toroidal surfaces corresponds to the value 
^center = 3.817. Interesting region containing both the cusps and the toroidal surfaces is plotted in detail at the left lower figure. 
Right lower figure shows the behaviour of the von Zeipel radius in the equatorial plane. 



characterizing rotating fluid configurations in terms of the LNRF orbital velocity is to introduce, so-called, von Zeipel 
radius defined by the relation 

{l-co£)p (15) 



LNRF 



which generalizes in another way as compared with [16] the Schwarzschildian definition of the gyration radius p. (For 
more details see [19].) 

In the case of marginally stable tori with £(r, 9) = const, the von Zeipel surfaces, i.e., the surfaces of £%(r, 9;a,£) = 
const, coincide with the equivelocity surfaces ^NR F (r, ^\ a A) = const. For the tori in the Kerr spacetimes, there is 



EVA(A-2fl£r) sin0 
A(A-a 2 sin 2 0)+4fl 2 r 2 sin 2 9 ' 



&{r,e;a,£) = (16) 
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FIGURE 2. Profiles of the equatorial orbital velocity of marginally stable tori in LNRF in terms of the radial Boyer-Lindquist 
coordinate (left column). For comparison, the profiles of the orbital velocity of Keplerian discs in Kerr spacetimes with the same 
rotational parameter a are added (right column). For thick discs, values of I = const are appropriately chosen; commonly, t = l m% 
is used giving the maximal value of the velocity difference in between the local extrema, and representing the limiting case of 
marginally stable thick discs. 



Topology of the von Zeipel surfaces can be directly determined by the behaviour of the von Zeipel radius (16) in the 
equatorial plane 

r(r 2 + a 2 )-2a(e-a) 



&(r, 9 = n/2;a,i) 



(17) 



The local minima of the function (17) determine loci of the cusps of the von Zeipel surfaces, while its local maximum 
(if it exists) determines a circle around which closed toroidally shaped von Zeipel surfaces are concentrated (see 
Fig. 1). Notice that the minima (maximum) of ^(r, 9 = n/2;a,£) correspond(s) to the maxima (minimum) of 



(9) 



LNRF 



r, 9 = 7r/2; a, £), therefore, the inner cusp is always physically irrelevant being located outside of the toroidal 
configuration of perfect fluid. Behaviour of the von Zeipel surfaces nearby the center and the inner edge of the thick 
tori orbiting Kerr black holes with a > a c ( t hick) = 0.99979, i.e., the existence of the von Zeipel surface with an outer 
cusp or the surfaces with toroidal topology, suggests possibility of strong instabilities in both the vertical and radial 
direction and a tendency for development of some vortices crossing the equatorial plane. We plan studies of these 
expected phenomena in future. 



5. DISCUSSION AND CONCLUSIONS 

It is useful to discuss both the qualitative and quantitative aspects of the phenomenon of the positive gradient of LNRF 
orbital velocity. In the Kerr spacetimes with a > a c (thick)> changes of sign of the gradient of ^W(r; a) must occur 
closely above the center of relevant toroidal discs, at radii corresponding to stable circular geodesies of the spacetime, 
where the radial and vertical epicyclic frequencies are also well defined. 
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FIGURE 3. (a) Positions of local extrema of I^krf (' n ^-L coordinates) for the marginally stable discs with £ = t ms in 

dependence on the rotational parameter a of the black hole, (b) Velocity difference = "f^tx — ^mfa as a f unct i° n °f tne 

rotational parameter a of the black hole for both the thin (Keplerian) disc and the marginally stable (non-Keplerian) disc with 



(b) 



In two interesting cases, behaviour of y^\r,a,l) is illustrated in Fig. 2; for comparison, profiles of the Keplerian 
velocity r; a) are included. With a growing in the region of a 6 (a c (thick) , 1). the difference AT 7 ^ = fmxl — 
grows as well as the difference of radii, Ar = r max — r m ; n , where the local extrema of y^\r;a,£) occur, see Fig. 3. 

In terms of the redefined rotational parameter (1 — a), the changes of sign of gradient of the LNRF orbital velocity of 
marginally stable thick discs occur for discs orbiting Kerr black holes with (1 — ci) < 1 — ^c(thick) ~ 

2.1 x 1(T 4 which 

is more than one order lower than the value 1 — fl c (thin) = 4-7 x 10~ 3 found by Aschenbach [4] for the changes of sign 
of the gradient of the orbital velocity in Keplerian, thin discs. Moreover, in the thick discs, the velocity difference, 
Af^ — Ymfx — "^mfj, is smaller but comparable with those in the thin discs (see Fig. 3). In fact, we can see that 
for a — > 1, the velocity difference in the thick discs A^^ k j w 0.02, while for the Keplerian discs it goes even up to 

A ^(t < hi!i) ~ °- 07 - 

are really huge velocity differences, being expressed in units of c. 
Following Aschenbach [4], we can define the typical frequency of the mechanism for excitation of oscillations by 
the maximum slope of the positive gradient of d"V^> /dr in the region of the changes of its sign by the relation 

dy/(<p) 

£2^, = 2k— 5— Ux. (18) 

The "oscillatory" frequency has to be determined numerically. We have done it for both Keplerian discs and the 
marginally stable discs with I = £ ms = const, see Fig. 4 and Table 1 . It is interesting that in Keplerian discs with 
a ~ 0.996, there is £2£ rit ~ £2r ~ £2y/3, i.e., £l r crit can be related to the resonant phenomena at the radius where 
the epicyclic frequencies are in the 3:1 ratio [4]. However, it is more correct to consider d'f^'/dR where R is 
the physically relevant (coordinate-independent) proper radial distance, as it is more convenient for estimation of 
physically realistic characteristic frequencies related to local physics in the disc. Then the critical frequency for 
possible excitation of oscillations is given by the relation 

d-fi.9) 



a* rit = 27T-^-| max , dR = Jg^dr = W -dr. (19) 

Of course, such a locally defined "oscillatory" frequency, confined to the orbiting LNRF-observers, should be further 
related to distant observers by the relation (taken at B-L coordinate r corresponding to [d'f^' / <3r) max ) 



&l = ^-{gtt + lag^ + aP-g^)^. (20) 

Similarly, an analogical relation can be written also for the critical frequency £2J, rit , giving the circular frequency 
£2^. Because the velocity gradient related to the proper distance is suppressed in comparison with those related 
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(c) (d) 

FIGURE 4. Critical "oscillatory"frequency for excitation of epicyclic oscillations, introduced by Aschenbach [4], as a function 
of the rotational parameter of the black hole in terms of both the B-L coordinate radius (£2£ rit ) and the proper radial distance (£2f rit )- 
(a) Keplerian discs, (b) Marginally stable (non-Keplerian) discs with constant specific angular momentum t = £ ms . (c) Comparison 
of critical frequencies for Keplerian fl^. it and non-Keplerian £2 "j? discs in terms of the proper radius, (d) Positive parts of the 
"coordinate" and "proper" radial gradient for a given value of the rotational parameter a. 

to the coordinate distance, there is i2^ rit < £2£ rit . The situation is illustrated in Fig. 4. Moreover, Fig. 4d shows 
mutual behaviour of the "coordinate" and "proper' ' radial gradient df^/dr and dY^/dR in region between 
the local minimum and the outer local maximum of the orbital velocity for an appropriately choosen value 
of the rotational parameter a. Characteristic frequencies / = D./2n, where Q. corresponds to the particular circular 
frequencies ^crit> ^°°> respectively, are given in Table 1. 

We can conclude that in constant specific angular momentum tori, the effect discovered by Aschenbach is elucidated 
by topology changes of the von Zeipel surfaces. In addition to one self-crossing von Zeipel surface existing for 
all values of the rotational parameter a, for a > a c (thick) the second self-crossing surface together with toroidal 
surfaces occur. Toroidal von Zeipel surfaces exist under the newly developing cusp, being centered around the circle 
corresponding to the minimum of the equatorial LNRF velocity profile. Further, the behaviour of von Zeipel surfaces 
in marginally stable tori orbiting Kerr black holes with a > a c ( t hi c k) strongly suggests a possibility of development of 
both the vertical and vortical instabilities because of the existence of the critical surface with a cusp, located above the 
center of the torus and the toroidal von Zeipel surfaces located under the cusp. The effect of "velocity gradient sign 
changes" can be very important as a trigger instability mechanism for oscillations observed in QPOs. Of course, further 
studies directed both to the theoretically well founded, detailed physical mechanisms for triggering of oscillations in 
the equilibrium tori with general specific angular momentum distribution, and the link to observations, are necessary 
and planned for the future. 

Finally, we would like to call attention to the fact that signs' changes of the radial gradient of orbital velocity in 
LNRF occur nearby the r — r^i orbit, while in the vicinity of the r = orbit, dY^^p/dr < for all values of a for 
both the Keplerian discs and the marginally stable toroidal discs with all allowed values of I. Clearly, the parametric 
resonance, which is the strongest one for ratios of the epicyclic frequencies ilv/£2p> = 3/2 works at the r = orbit, 
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TABLE 1. Characteristic frequencies in units (M/M Cv )~ l Hz (M/M@ is the mass of 
the Kerr black hole in units of mass of the Sun.) related to the circular frequencies 
£2 = 2nf defined in the text are given for appropriate values of the black hole spin. 
Maximal values of the frequencies related to the stationary observer at infinity are bold- 
faced. 
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while its effect is much smaller at the radius r = r^,- 1 with i2y / £2r = 3/1 [9] . Therefore, the forced resonance, triggered 
by the changes of d^Kgp/dr, will be important for the 3:1 resonance. Notice that the forced resonance at r = r$-\ 
can generally result in observed QPOs frequencies with 3:2 ratio due to the beat frequencies allowed for the forced 
resonance [20]; however it seems to be irrelevant in the case of microquasars, as all observed frequencies lead to the 
values of the rotational parameter a < fl c ( t hick) as shown in [13]. 
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